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Abstract. We present integral field spectroscopy of the LBV candidate stars B 416 and v 532 in the local group galaxy M 33. 
B 416 is surrounded by an elongated ring-like nebula, which has a projected radius of 20 x 30 pc. From the datacube we create 
ionization and radial velocity maps of the nebula. The excitation of the gas decreases towards the outer part of the ring, while 
the inner part of the nebula is filled with a more excited gas. In the EW direction the ring is seen to expand with a maximum 
projected velocity amplitude of about 40km/s. The eastern part approaches the observer. We estimate the nebula dynamical 
lifetime ~ 8 • 10^ years. It could be a residual MS bubble, which indicates a main-sequence or pre-LBV status of the star. 
We classify B 416 as an “iron star” or B[e]-supergiant. In v532 an elongated nebula has been marginally detected. The total 
projected size of the nebula along the main axis is 30 pc, and the total radial velocity gradient is 44± 11 km/s. v 532 exhibits both 
strong photometric and spectral variability. At the time of our observations it was in an intermediate brightness state with a rich 
nitrogen spectrum. We classify v 532 as an LBV, showing LBV Ofpe/WN transitions. We stress the importance of integral 
field spectroscopy as the optimal technique for studying nebulae and the evolution of LB V-like stars in nearby galaxies. 
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Introduction 

Analysis of the brightest supergiants as individual stars in ex¬ 
ternal galaxies provides a unique tool for determining the prop¬ 
erties of young stellar populations in their host galaxy. In par¬ 
ticular, nearby galaxies are ideally suited to studying the most 
massive stars systematically, both during their normal evo¬ 
lution in different environments, as well as during the short 
and unstable stages of evolution as OB supergiants, rare hy¬ 
pergiants, Luminous Blue Variables (LBVs), Ofpe/WN stars, 
B[e]-supergiants (Kudritzki 1998; Lamers et al. 2001), or dur¬ 
ing their latest stages of evolution as peculiar massive binary 
objects with relativistic jets and accretion disks like SS 433 
(Fabrika and Sholukhova 1995). All of these stars are young 
and massive, and they are among the brightest stars in their 
host galaxy. They share some main spectral properties like a 
blue continuum and strong emission lines, and they often show 
a surrounding nebula. In the Milky Way a considerable fraction 
of such stars may be hidden by dust extinction in the galactic 
plane, which is why they are observed best in nearby galaxies. 

The evolutionary connections between the various classes 
of massive stars, such as the OB-stars, B[e]-supergiants, LBVs, 
Ofpe/WNL, and Wolf-Rayet stars are not well understood. 
Systematic surveys of LB V-like candidates in massive galaxies 
(Calzetti et al. 1995, Massey et al. 1996, Corral 1996, Fabrika 
and Sholukhova 1999) provide a basis for isolating these ob¬ 
jects for subsequent follow-up and quantative spectroscopy 


(Massey et al. 1996, Sholukhova et al. 1997, 1999) to reveal 
the basic properties of these objects. In M 33 there are 4 con¬ 
firmed LBV stars and about a dozen LBV candidates, but only 
one known Ofpe/WN star, and about 5 new candidates for this 
latter class of objects (Massey et al. 1996). Not a single SS433- 
type star has been found in M 33, and also no B[e]-supergiant 
has been confirmed to date. 

LBVs are extremely massive stars (M > 40 Mq) at 
an unstable stage of evolution with a mass loss rate of ~ 
10“^ MqIj. They show strong spectral and photometrical vari¬ 
ability on time-scales from months to years (Humphreys and 
Davidson 1994; Lamers et al. 1998a). LBVs represent an insta¬ 
bility stage after the main sequence. It is not clear whether there 
is a relation between LBVs and B[e]-supergiants (Conti 1976), 
since their spectra often look practically the same, but the lat¬ 
ter stars do not show strong variability (Zickgraf et al. 1986; 
Lamers et al. 1998b). Spectral properties of B[e]-supergiants 
can be interpreted by a non-isotropical mass-loss occurring in 
the equatorial plane (Zickgraf et al. 1985; Zickgraf et al. 1986). 
They are probably a subclass of rapidly rotating post main- 
sequence stars. One may expect an evolutionary connection be¬ 
tween the B[e]-supergiants and the LBVs. 

Important information can be derived from morphological 
studies of the nebulae surrounding such stars. Because of their 
heavy mass loss, LBV stars are expected to show extended en¬ 
velopes, which consist of ejected matter and swept-up interstel¬ 
lar material. The same is true for B[e] supergiant stars. Direct 
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imaging and kinematic studies help to shed light on the mass 
loss history and on the physics of the stellar wind. 

Numerical models of the nebular expansion around massive 
stars (Garcia-Segura et al., 1996) have shown that those nebu¬ 
lae reach diameters of up to 20 - 40 pc at the main-sequence 
and pre-LBV stages. Such scales correspond to ~ 10" at the 
distance of M 33 (3"5/pc). LBV nebulae are smaller because of 
a shorter evolution time. Practically all known LBV stars have 
circumstellar shells (Humphreys and Davidson 1994). Typical 
LBV nebula sizes are in the range of 0.1 - 4 pc, expansion ve¬ 
locities 15 - 100 km/s, and dynamical times are in the range 
100 - 5 • 10^ years (Nota 1999; Figer et al. 1999; Weis 2003). 
There are possibly larger gas nebulae around LBV stars, for ex¬ 
ample a ring-like nebula around S Dor has a size up to 70 pc 
(Weis 2003). 

Spectroscopic analysis of the star is often hampered by 
the presence of strong nebular emission lines which contami¬ 
nate the stellar spectrum, and is complicated in crowded fields. 
For our targets we choose the technique of Integral Field 
Spectroscopy (IFS, sometimes also called 3D Spectroscopy) in 
order to optimize the separation of the stars from their crowded 
environments. It provides both improved subtraction of nebu¬ 
lar contamination from the stellar spectrum and spatially re¬ 
solved nebular spectra. One can make extensive use of the full 
2-dimensional information contained in 3D data, which can be 
also described as stacks of very narrow bandwidth filter images 
(datacubes). 

Here we report on results of IFS observations of two LBV 
candidates in the Local Group galaxy M 33, the stars B 416 and 
v532. The first star was designated as blue star N416 in the 
photographic survey by Humphreys and Sandage (1980). It was 
also identified as an LBV candidate in follow-up spectroscopy 
of the LFV-brightest stars of M33 (Massey et al. 1996). In a 
program searching for SS 433-like stars in M 33, Shemmer and 
Leibowitz (1998) detected periodic microvariations in the pho¬ 
tometry of B416. These oscillations were confirmed as peri¬ 
odic and coherent by Shemmer et al. (2000), who classified 
the star as LBV based on the spectral similarity in comparison 
with typical LBV spectra at a quiescence stage. The coherent 
brightness oscillations may become quite important for our un¬ 
derstanding of the physics of LBV stars; however, B 416 does 
not show the large amplitude variability typical for LBVs. 

The second object of this paper is the variable blue star 
V 532 (Romano 1978; Artyukhina et al. 1995), which is known 
to exhibit the characteristic variability of LBVs, i.e. an ampli¬ 
tude of « 1™ on a time-scale of years (Kurtev et al. 2001). 
Short-scale and quasi-periodical variabilities have also been 
found in this star (Sholukhova et al., 2002). Strong emission 
lines, spectral and photometrical variability (Szeifert 1996; 
Fabrika 2000; Sholukhova et al. 2002; Polcaro et al. 2003) in¬ 
dicate that V 532 is an LBV object. 

Observations and data reduction 

Observations of v532 and B416 were carried out with the 
Multi-Pupil Fiber Spectrograph MPFS' (Afanasiev 1998) at the 

* http://www.sao.ru/''gafan/devices/mpfs/mpfs_main.htm 


prime focus of the 6m telescope BTA (Russia). The integral 
field unit (IFU) of 16x 15 square spatial elements (’’spaxels”) 
covers a rectangular region of 16" x 15" on the sky (1" per el¬ 
ement). Optical fibers transmit the light from the 240 spaxels 
of the IFU together with 16 additional fibers located at a radius 
of 4'.5 outside of the IFU for the purpose of sampling the sky 
background far away from the field. In nearby galaxies like in 
M33 the sky fibers will still be located within the galaxy, thus 
we did not use the sky fibers for the background subtraction. 
The entire set of 256 fibers is reformatted to form a pseudo-slit, 
whose emerging light is dispersed and projected onto the focal 
plane of the spectrograph. We used a SITe TK1024 backside- 
illuminated CCD with 1024x1024 pixels and a pixel size of 
2Apm. 

Two 20-minute exposures of v532 were taken on 
September 18, 1998 with a reciprocal dispersion of 1.3 A/pixel 
(FWHM Si 3.5 A) in the spectral range 4470 - 5800 AA 
with a seeing of 1.5-2 arcsec FWHM. B 416 was observed on 
September 28, 1998 in two 20-minute exposures with a recip¬ 
rocal dispersion of 2.6 A/pixel (FWHM ss 7 A) in the spectral 
range of 4250 - 6900 AA with a seeing of 3 arcsec FWHM. 

3D observations of B416 were also carried out with 
INTEGRAL^ at the 4.2m WHT, La Palma (Arribas et al. 1998) 
on January 18, 2001. This instrument is a bare fiber bundle 
type of 3D spectrograph. Contrary to MPFS, the IFU has no 
lensarray in front of the fibers, resulting in a hexagonal pack¬ 
age with small gaps between the circular fiber apertures. We 
used the SB 1 bundle, which has 205 fibers and covers a rect¬ 
angular FOV of 7.8x6.4 arcsec^ with a projected fiber diame¬ 
ter of 0'.'45. The fiber bundle is coupled to the bench-mounted 
fiber spectrograph WYFFOS on the Nasmyth platform of the 
WHT. The detector is a thinned, backside-illuminated IKxlK, 
24 pm pixel CCD. We used an R300B grating, giving a spec¬ 
tral coverage of 5800 A and a spectral resolution of 11 A at 
a reciprocal dispersion of 5.8 A/pixel and took 2 exposures of 
1800 sec each. The conditions were clear with a seeing of 0.9- 
1.0" FWHM. 

Data reduction for both the MPFS and INTEGRAL data 
sets was performed with the P3d package (Becker 2002), which 
was originally developed for the PMAS instrument (Roth et al. 
2000), but has also been shown to work well for other instru¬ 
ments. The raw CCD frames were bias and dark subtracted, 
and cleaned from cosmic ray events. The spectra extraction 
was performed using a profile fitting algorithm, which simul¬ 
taneously solves for cross-talk between adjacent spectra and 
straylight patterns on the detector. P3d incorporates a robust 
algorithm to reliably determine the geometry of the spectra. 
After extraction a wavelength calibration was applied using arc 
exposures. The wavelength-dependent fiber-to-fiber transmis¬ 
sion variation was calibrated by means of twilight Hatfield ex¬ 
posures. At this final stage of data reduction a 2-dimensional 
image of stacked spectra appears. Monochromatic maps or an 
entire datacube may be subsequently derived from this stacked 
spectra format. The final analysis of the fully reduced 3D 
data, i.e. generation of monochromatic maps at selected wave¬ 
lengths, co-adding flux within a digital aperture, sky subtrac- 

^ http://www.ing.iac.es/~bgarcia/integral/html/integraLhome.html 
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Fig. 1. (a) MPFS map of B416 at 6000 A (spatial sampling 1 arcsec, the FOV is 16x15 arcsec^) and the same field (b) from the 
HST WFC2 image. An enlarged fraction of this field from the INTEGRAL in the same wavelength (c) and the corresponding 
enlarged HST field (d). Note the different orientation of the MPFS and INTEGRAL fields. 


tion, etc., was performed with the MONOLOOK tool (Becker 

2002). 

The full 2-dimensional spatial information allows one to 
accurately model the background of resolved or unresolved 
stars and gaseous emission at any wavelength. A detailed de¬ 
scription of this technique is given in Roth et al. (2003). 
Application of image deconvolution techniques is able to dis¬ 
entangle the different components from severely crowded stel¬ 
lar fields with significant contamination from nebular emission. 
We used the cplucy two-channel algorithm, which is available 
in IRAL. A more detailed account of this technique is given in 
Becker (2002) and Becker et al. (2003). 

In Lig. 1 we present continuum maps of B 416 from a wave¬ 
length region near 6000 A of MPLS and INTEGRAL, respec¬ 
tively. Lor each ILU map the same LOV as observed with HST^ 
is shown for comparison. While the MPLS map (a) represents 
the intensity distribution exactly over the square spaxels of the 
lens array, the INTEGRAL map (c) was interpolated from the 

^ Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope 
Science Institute. STScI is operated by the Association of Universities 
for Research in Astronomy, Inc. under NASA contract NAS 5-26555. 


205 data points of the fiber bundle (overlay). Note the patches 
of low intensity near the upper left corner in (c) indicating the 
presence of dead fibers in the bundle. 

While the MPLS observations suffered from rather poor 
seeing, the conditions with INTEGRAL were more favourable 
and yielded much better spatial resolution. In Lig. 1(a) the im¬ 
age of B 416 is blended with the star cluster located from 2"5 
to 5"5 to the SW. In Lig. 1(c) the cluster is clearly separated 
from B 416 and marginally resolved into stars. Three other iso¬ 
lated stars seen in the HST frame to the W of B 416 are also 
marginally present in the UNTEGRAL map. 

Results 

The nebula around B416 

Ligure 2 presents monochromatic images of B 416 in the emis¬ 
sion lines of Ha, li/3, He I 45876, [SII] 46717 H- 46731, 
[Nil] 46583, and [GUI] 45007 from the MPLS observations. 
The images were produced in +2.5 A bands from the cen¬ 
tral pixels of the lines (7.5 A bands) as relative fluxes. 
Corresponding continuum images were substracted. There is 
a shift in the monochomatic pictures depending on wavelength 
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Fig. 2. Monochromatic MPFS maps of B416 taken in central 7.5 A bands of different emission lines, continuum subtracted. 
Intervals of relative fluxes (counts) used in the maps are shown. The position of the star is indicated by a white cross. A prominent 
arc-like nebular feature is visible to the left. The star itself shows emission in the Balmer lines and He I. 


due to atmospheric dispersion. Its value is O'.'42/IOOOA (Cf.'hS 
between Ha and [OIII] 45007) towards the NW with decreas¬ 
ing wavelength. The central star position is indicated in the fig¬ 
ure by crosses in each image taking the atmospheric shifts into 
account. 

A ring-like nebula is clearly seen in the emission line maps 
of Fig. 2. It is non-symmetrical and presents a different mor¬ 
phology in lines of different excitation. The star itself is a 
source of emission in permitted lines. The stellar emission in 
Ha is particularly strong and gives rise to a high level of con¬ 
trast, which is why the nebula appears faint at this wavelength. 
There is no apparent forbidden line emission coincident with 
the PSF of the central star. The inner parts of the nebula are 
filled with faint and diffuse emission in [OIII] 44959,5007. 

In order to disentangle nebula and stellar PSF components 
in the spectral domain we attempted a gaussian decomposi¬ 
tion of the nebular and stellar emission line profiles in Ha and 
HjS. Since the nebular line is intrinsically narrow, but the stellar 
emission line profile is expected to be broader under the influ¬ 


ence of a strong stellar wind, this approach seemed to provide 
a useful criterion. 

As a best guess we estimated the width of the nebular 
emission in Ha and [Nil] 46548,6584 in regions far away 
from the star, yielding a gaussian FWHM ss 6.8 A, which is 
in fact identical to the spectral resolution. We then fitted the 
observed spectrum near Ha by three gaussians for Ha and 
[Nil] 46548,6584, whose widths are equal to that of the instru¬ 
mental profile, and by a fourth component for the stellar Ha 
emission, whose FWHM was determined from the fit. In fact, 
for all spaxels close to the centroid of the star, this latter compo¬ 
nent makes a non-negligible contribution to the spectrum and 
has an FWHM 12.2 A. The broad component contributes to 
55 % of the total Ha emission inside of the centroid of the star, 
but it is absent in spaxels far away from the star. The observed 
spectrum in a spaxel close to the star, the three narrow and one 
broad gaussian fits, and the sum of the fits are plotted in Fig. 3. 
This spectrum was extracted in a spaxel not far from the star to 
show both broad and narrow components of Ha. The flux units 
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are the same as in Fig. 2, whose images were produced in 7.5 A 
bands. Similarily, the decomposition of Hj6 into broad and nar¬ 
row profiles resulted in a broad line width of FWHM ss 10 A. 
This value agrees well with the value found in Ha and with the 
corresponding intrinsic velocities. 



Fig. 3. An example of the Ha and [Nil] /16548,6584 lines 
(solid line) in a region close to the star. They are fitted with 
four gaussian profiles (dashed lines) — broad and narrow Ha 
lines and two narrow [Nil] lines. A sum of the gaussian profiles 
is shown by a long-dashed line 

Figure 4 demonstrates a spatial representation of the spec¬ 
tral decomposition into narrow and broad profiles, where the 
broad component is associated with the star and is perfectly 
coincident with the stellar PSF. In the nebular analysis, which 
follows below, we only use the narrow B aimer line component. 
These results suggest that the broad Balmer emission lines are 
indeed most likely formed in the atmosphere of the star. We 
must stress both that an expanding nebula near the star (with 
an angular extent below our spatial resolution) cannot entirely 
be ruled out, and that even the stellar atmosphere can contribute 
to the narrow emission. 

The morphology of the ring-like nebula presents an almost 
perfect elliptical ring (Fig. 2), which is however interrupted on 
its NW side by a distinct gap. The overall surface brightness 
distributions in Ha, [Nil], and [SII] are fairly similar, but show 
differences in the exact size and morphological details. The sur¬ 
face brightness distribution in [OIII] nebula is strikingly differ¬ 
ent from those in the Ha, [SII], and [Nil] lines. The extended 
faint emission in the N, W parts and very bright SW regions 
are dominated by high excitation gas contrary to those parts of 
opposite directions. The H/l nebula pattern (and Hy, which is 
not shown here) is somewhere intermediate between those of 
high excitation ([OIII]) and low excitation ([SII] and [Nil]). 

To illustrate the whole nebula around B 416, we present in 
Fig. 5 a small (roughly 1 arcmin) subfield from Ha and [OIII] 
images of M33 taken by Massey et al. (2001) in the course of an 
imaging survey of Local Group galaxies, though the Ha neb¬ 
ula was first detected by Shemmer et al. (2000). The images of 
Massey et al. (2001) were obtained with the CTIO and KPNO 
4-m telescopes using the Mosaic CCD cameras. Narrow-band 


filters were used with bandwidths of 55 A in [OIII] and 80 A in 
Ha. The MPFS field and corresponding emission line isophotes 
of the narrow Ha line and [OIII] 45007 are overplotted. Each 
step in the isophotes corresponds to 200 counts, the same rel¬ 
ative units as used in Fig. 3. The maximum isophote level is 
1900 counts in Ha and 900 counts in [OIII] 45007. The MPFS 
maps are in excellent agreement with these direct images. 

The [OIII] emission around B 416 is faint. It is concentrated 
on the SW side of the nebula, where an extended [OIII] region 
does exist, which is only partly covered by the MPFS field. 
The most distant SW part of the nebula (just outside of the 
MPFS map) has about the same shape in [OIII] and Ha (Fig. 5). 
This fact, together with the specific shape of the [OIII] nebula, 
may point towards an explanation that the most distant S W part 
of the whole nebula is shaped not only under the influence of 
B 416, but also by a cluster of hot stars. Two compact clusters 
of stars (Fig. lb) are located in these regions, coinciding with 
the [OIII] bright emission of the nebula, where, conversely, 
[SII] and [Nil] are relatively faint. Yet another high excitation 
knot is seen within the diffuse [OIII] region, just 3" SW of the 
edge of the MPFS FOV (Fig. 5). One could speculate that the 
bright [OIII] emission is due to photoionization from hot stars 
in these clusters, whose relation, however, to the ring-like neb¬ 
ula is not immediately clear. 

A similar situation is observed in the nebula DEM L106 
around EMC B[e]-supergiant Hen S22 (Chu et al. 2003). The 
HII region N30B with a cluster of hot stars is located inside the 
large nebula (~ 30 pc in radius). A bow-shock-like halo around 
N30B is oriented towards Hen S22. It is formed as a result of 
interaction with the stellar wind of Hen S22. The halo reflects 
the supergiant radiation. One could conclude that the most dis¬ 
tant SW region of the ring-like nebula of B 416 is probably 
affected in a similar fashion in the environment of the hot star 
cluster. 

Eor a qualitative analysis of the excitation structure of the 
nebula we present “ionization” maps of B416 in Eig. 6. Low 
excitation regions are traced by the [SII]/Ha map, whereas high 
excitation regions appear in the [OIII]/H/3 ratio map. In both 
cases we used the Balmer line maps derived from the narrow 
line components as described above. One can see from a com¬ 
parison of Eig. 6, 2, and 4 that over a wide range of position 
angles, the inner parts of the ring-like nebula are dominated 
by high excitation, as opposed to the outer rim of the ring, 
which is most pronounced in the low ionization radiation of 
[SII]. This observation is qualitatively in line with the expecta¬ 
tion of finding (i) photoionization predominantly near the star 
and (ii) shock excitation from the interaction of the expanding 
shell with the interstellar medium at the edge of the nebula. 

We measured the diagnostic line ratios of Ha/[NII], 
Ha/[SII] (sum the doublets) and I(6717)/I(6731) for [SII], and 
plotted the results in the diagnostic diagrams of Sabbadin & 
D’Odorico (1976). Our mean values for the nebula around 
B416 are as follows: Ha/[NII]= 3.24 + 0.39, Ha/[SII]= 2.21 + 
0.25, and I(6717)/I(6731)= 1.61 +0.11, where the r.m.s. values 
are for individual resolution elements. The mean line ratios are 
quite typical for HII regions; however, the ratios are not con¬ 
stant accross the nebula. It is directly seen in Eig. 2 that the rim 
is more distant in [SII] than in [Nil] and Ha. The outer parts of 
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Fig. 4. Maps of B416; narrow H/3 component (left panel), narrow Ha (middle), and broad Ha (right). The stellar centroid is 
indicated by a cross in each frame. 
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Fig. 5. Ha (left) and [OIII] /15007 (right) images of the B416 region taken by Massey et al. (2001) with the MPFS field and the 
line izophotes imposed on the images. 


the ring, at least in the NES directions, consist of a gas of lower 
excitation. This conclusion is confirmed by the [SII]/Ha ratio 
map (Fig. 6). The ratio is greater than 0.4 in the outer parts 
of the ring, which indicates a collisional excitation of gas as 
was confirmed by Smith et al. (1993) for SNRs in M33. The 
electron density (Osterbrock 1989) over the whole nebula was 
found to be at the low density limit, n^ < 100 cm“^, with an 
average value of «« ~ 10 cm“^. 

The lower left panel of Fig. 6 shows a map of the Ha/H;6 
ratio, whose canonical value based on recombination theory is 
expected to be si3. A systematic increase of this value is seen 
along the E rim of the ring-like nebula. Most strikingly, how¬ 
ever, the ratio increases to a value of 5 at the N border of the 
map. This peculiar region happens to coincide with the afore¬ 
mentioned gap in the structure of the ring (Fig. 5), which is thus 
explained as extinction from dust. 

A point-by-point analysis of radial velocities in Ha, 
[OIII] 45007, and [SII]/16117 lines was performed by fitting 
gaussians to each line. The internal accuracy of the instrumen¬ 
tal wavelength calibration was increased by measuring differ¬ 
entially against the [OI] 45577 and 46300 night sky lines. The 
resulting accuracy of the radial velocity measurements is about 
10 km/s per spatial element (1" ) in those places where the 


nebular emission provides sufficient intensity (Fig.2). The fi¬ 
nal radial velocity maps were averaged over 3 pixels, and the 
resulting radial velocity accuracy is better than 10 km/s per res¬ 
olution element. 

The lower right panel of Fig. 6 presents a radial velocity 
map as derived from the deblended [SII] 46717 emission line. 
The overplotted isophotes were derived from the total intensity 
of the doublet [SII] 446717,6731 (cf. Fig. 2). We chose this 
emission line because it is bright and ubiquitous over the en¬ 
tire FOV, thus minimizing the uncertainty of our measurement. 
The overall velocity field pattern, however, is the same for all 
of the bright nebular lines, and shows a systematic variation 
with a peak-to-peak amplitude of about 30 - 40 km/s. There 
are two conspicuous features of extreme radial velocities: an 
area of roughly 9" x5" across the E rim of the nebula has a 
minimal velocity of -226 km/s, whereas an area of the same 
extension and at the same distance to the W of the star has a 
significantly more positive radial velocity of ss -184 km/s. In 
the N-S direction perpendicular to these two patches the radial 
velocity is constant (AV^ < 5 km/s) and zero with respect to the 
average, except for the S corner of the field, where the velocity 
is about -hIO km/s more positive than the average. 
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Fig. 6. Ionization maps of B416 obtained as ratios of line intensities of [OIII] /15007/HjS(narrow), [SII] (/16717 + 
/l6731)/HQ;(narrow), a map of Ha/H/6 ratio and radial velocity map obtained in [SII] A61\l emission with relative flux izophotes 
of the sum of the lines of the [SII] doublet (compare with Fig. 2). Crosses indicate 4 selected regions where mean radial velocities 
were measured (see text). 


We derived mean radial velocities of Ha, [SII] 46717, and 
[OIII] 45007 lines in four selected E, N, W, and S regions 
marked by crosses in Fig. 6. The mean velocities of the Ha, 
[SII], and [OIII] lines respectively are in the E-region -239 ± 
3, -226 ± 2, -246 ± 8 km/s; in N-region -226 + 2, -210 ± 
3, -232 + 5 km/s; in W-region -212 ± 3, -184 + 4, -216 ± 
3 km/s, and in S-region -220 ± 3, -197 + 7, -213 ± 4 km/s. 
The errors of the mean values depend both on brightness and 
homogeneity of a line emission in these selected regions. Some 
systematic differences in radial velocities between these three 
lines may appear because they were formed under different 
physical conditions (unresolved gas clouds and structures). The 


measurements in the selected regions confirm the total radial 
velocity gradient observed in the nebula in the E-W direction. 
It is /S.Vr(EW) = 33 ± 8 km/s over these three lines, and is 
AVriEW) = 34 ± 4 km/s in Ha + [SII] 46717, which are the 
most similar both in physical conditions and in nebula patterns 
(Fig. 2). 

The star’s radial velocity in Ha (the total profile) is 
-231.1 ±3.5 km/s. It is intermediate between the E and W 
opposite features of extreme radial velocities in Ha and closer 
to radial velocity of the E region in accordance with location 
of the star inside the ring-like nebula (Fig. 2). The stellar ra¬ 
dial velocity may be considered as the systemic velocity of the 
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whole complex (star + nebula). However the total stellar Ha 
line profile may be distorted because of the contribution (50 %) 
by the broad stellar component and stellar wind. 

We summarize our findings that the morphological asym¬ 
metry of the nebula in the E-W direction is accompanied by 
an asymmetry in the radial velocity pattern; the E part of the 
ring approaches (» -15 km/s), while the W part receedes 
(-1-15 -1- h- 20 km/s). There is no apparent asymmetry in the N-S 
direction. The asymmetry of the nebula in the E-W direction is 
accompanied also by an asymmetry in the gas excitation, while 
the approaching E part of the ring is dominated by low exci¬ 
tation emission (a collisional excitation in the rim), but the re- 
ceeding W part is dominated by high excitation gas (a radiative 
excitation by the central star). 

A nebula around v532 (Romano’s star) 

Our MPES images of v532, another LBV-candidate in M33, 
were taken merely in the spectral range 4470 - 5800 AA (spec¬ 
tral resolution » 3.5 A). In the absence of Ha and [SII] we had 
to resort to an analysis in H/l, since no [OIII] 44959, 5007 were 
detected in the spectrum. Although at first glance the H/l and 
continuum maps of this object appeared point-like, a suspicious 
slight elongation of the stellar image in H/3 with P.A. ~ 35° led 
us to perform a gaussian decomposition of the spectral profile 
into a broad and a narrow component along the lines of the ex¬ 
ercise with B 416 as described above. As before, we stress that 
besides an extended nebula, a narrow component may either 
originate in the stellar atmosphere or in a nearby unresolved 
nebula, or in both. We obtained a good fit from two gaussians 
with PWHM(narrow) = 4.1 A and PWHM(broad) = 11.1 A. 

The resulting maps of v532 in broad H/i, narrow H/l, and 
in the nearby continuum are shown in Pig. 7. Note that the 
narrow-line component is a factor of 4 brighter than the broad 
component. In order to assess the possible spatial extent of 
these surface brightness distributions, which are close to un¬ 
resolved point sources, we performed a Moffat fit to the three 
maps and inspected their radial intensity profiles for compar¬ 
ison. Pigure 8 shows that the narrow component is more ex¬ 
tended (371 PWHM) than the broad component (276 PWHM), 
which in turn has the same spatial extent as the continuum 
intensity profile. The difference between the broad and nar¬ 
row component intensity profiles becomes significant at a scale 
greater than the seeing value, ss 2" PWHM. The error bars in 
the intensity profiles are about the same. In five consequent 
bins from 3" to 5" the difference between the broad and nar¬ 
row component energy distributions varies from 3.5 to 2.5 stan¬ 
dard deviations. We conclude that an extended nebula has been 
marginally detected in the narrow emission line component of 
the object. 

A radial velocity analysis of the narrow Hj6 component sup¬ 
ports this view. We measured the radial velocity in each spaxel 
differentially against the night sky [OI] 45577 emission line as 
a reference. The central part including the star has a velocity 
« -220 km/s. The SW part of the nebula shows on average 
positive relative velocities, while the velocities of the NE part 
are mainly negative. In order to check this finding against any 



Fig. 8. Radial intensity profiles of v532 maps in narrow (filled 
symbols) and broad components (open circles) of the H/i line 
profile. 
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Fig. 9. V 532 nebula radial velocity map derived in Hj6 line. The 
cross marks the position of the star. 

potential errors from the spectral decomposition, we repeated 
the velocity measurements by fitting single gaussian profiles 
to the spectra which had been coadded from several spaxels in 
three representative regions (Fig. 9): center (12 spaxels), NE (9 
spaxels), and SW (4 spaxels). This procedure increases the S/N 
ratio in the resulting compound spaxels. The derived H/? radial 
velocities are -219±2.7km/s in the central part, -191 + lOkm/s 
in the SW, and -235 ±11 km/s in the NE part. 

The analysis confirms the presence of a radial velocity gra¬ 
dient across the nebula; the SW side of the nebula recedes 
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Fig. 7. V 532 maps in broad (left), narrow H/6 (middle), and in the nearby continuum (right panel). Crosses mark the centroid 
of the star. 


and the NE side approaches. The total radial velocity gradient 
amounts to 44 ± 11 km/s. The total angular size of the nebula is 
9" (ss 30 pc in projection) in the NE-SW direction. Erom the 
heterogeneous velocity distribution in Eig. 9 one may suspect 
that the system is more complex than just the simple picture of 
a bipolar nebula with receding and approaching lobes. Given 
the limited spatial resolution of our data we therefor forgo any 
further analysis. 

The Spectra of B 416 and v532 

3D and longslit spectra of B 416 and v 532 are shown in Eig. 10. 
Three rectified spectra of B 416 from different instruments are 
plotted with an offset for clarity. The top spectrum was ob¬ 
tained with the BTA long-slit spectrograph on January 19, 2001 
(a spectral resolution EWHM » 7 A). The middle spectrum 
was obtained by coadding spaxels with a digital aperture of the 
MPES datacube, observed on September 28, 1998 (EWHM « 
7 A). The bottom spectrum comes from the INTEGRAL dat¬ 
acube, which was observed on January 18, 2001 (EWHM « 
11 A). The lower panel in Eig 10 presents two v532 spectra: 
the top spectrum was obtained with MPES on September 18, 
1998 (EWHM » 3.5 A) and the lower spectrum from the long- 
slit spectrograph, observed on July 12, 1999 (EWHM » 3.5 
A). 

A comparison of the long-slit (LS) B416 spectrum with 
its 3D spectroscopy counterparts (MPES, INTEGRAL) demon¬ 
strates, that the former is severely affected by nebular contami¬ 
nation, e.g. [GUI] 44959,5007. The strongest line [OIII] 45007 
is superposed on a blend of Ee II emission lines in B 416 spec¬ 
trum. The nebular lines, however, are absent in the MPES and 
INTEGRAL plots, because the process of 2-dimensional back¬ 
ground correction accurately subtracts the local contribution of 
this light. In general, the flux distribution along a slit is not re¬ 
ally adequate to estimate the background at the location of the 
object. 

We inspected the spectra for spectral variability. According 
to Shemmer et al. (2000), B416 shows only weak photomet¬ 
ric and spectral variability, which is more or less in line with 
our observations taken essentially at two epochs in 1998 and 
2001 (the LS and INTEGRAL spectra were coincidently taken 
only one night apart from each other). Considering the differ¬ 


ence in spectral resolution, the LS spectra are indistinguishable 
from INTEGRAL, despite the systematic [OIII] contamination 
at 444959,5007. There are some small differences between LS 
and MPES spectra. The situation is quite different in the case 
of V 532, which shows a pronounced onset of P Cyp profiles in 
He I between 1998 (MPES) and 1999 (LS). 

The most prominent lines in both stars are the strong emis¬ 
sion hydrogen lines with broad wings and narrower He I lines. 
He I 44922,5015 lines in B416 are blended with the brighter 
Eell 44923,5018 lines, but the single He I 45876,6678 lines 
are clearly present in the red part. 

The main difference in the spectra of these stars consists, 
however, in the chemical composition of their atmospheres. 
There are many Eell lines in B416 and Nil lines in v532. In 
B 416 besides Ee II, Ee III, and He I lines, the following ions are 
also detected: Si II, Cr II, NII, Ti II, SII, CII, and also neutral 
elements such as Til and Mg I are present in the atmosphere. 
The Ee II wavelengths are indicated in the B 416 spectrum (un¬ 
labeled in Eig. 10), while the same is true for some NII lines in 
the v532 spectrum. 

The star v532 exhibits appreciable spectral and photomet- 
rical variability (Eabrika 2000; Kurtev et al. 2001; Polcaro et 
al. 2003). In some epochs strong PCyg components appear 
in He I or the hydrogen lines. Sometimes the He I emission 
component disappears, sometimes it becomes very bright, and 
the He II 4 4686 line appears in emission. In our MPES spec¬ 
trum the He I lines are bright with weak blue-shifted absorp¬ 
tion components. A PCyg line profile is more clearly seen in 
He I 44922, split by 220 km/s. A faint He II 44686 line is pos¬ 
sibly present (Pig. 10) with P Cyg-like profile. The spectrum is 
strongly blended by weak emission lines of NII, N III, Ee II, 
Ee III, CII, C III, Si II, Si III. Paint emissions of Ti, Cr, and S 
are also observed. The most numerical are nitrogen lines. A 
typical NII broad emission blend ss 5000 A is observed. The 
signal-to-noise value in the MPES spectrum is at least 30 in 
4700 A and 5300 A regions, and is no less than 40 around 
5000 A. 

Both stars possess extended atmospheres, as is apparent 
from the presence of the forbidden [Pelll] and [Pell] lines, the 
brightest of which are marked in the spectrum. We have also 
detected the [Nil] 45754.8 forbidden line in both stars (not 
shown in Pig. 10). 
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Fig. 10. Normalized spectra of B416 (top panel), and v532 (bottom panel). B416 spectra are from INTEGRAL (bottom), 
MPFS (middle), and the LS-spectrograph (top), v 532 spectra are from the LS-spectrograph (bottom) and MPFS (top). Numerous 
Fell lines in the B 416 spectrum and Nil lines in v 532 spectrum are not labeled in order not to confuse the plot. 


Discussion and conclusions 

One may suggest that B416 is surrounded by an expanding 
ring-like nebula, whose overall morphological, kinematic, and 
gas excitation properties are compatible with an outflow from 
the star, which happened to be confined to a plane. If the ring¬ 
like appearance were due to an increased column-density near 
the edge of an assumed spheroidal nebula, we would not expect 
to find the two patches of extreme radial velocities confined to 
the approaching and receding parts of the ring. It is possible 
that the plane of the torus is coincident with the equatorial plane 
of the star. The interior region of the nebula is filled with hot, 
high excitation gas. The pronounced low excitation emission 
near the edge of the ring indicates a shock ionization of material 
running into the interstellar medium. 

Fitting an ellipse to the entire Hq, images of Massey et al. 
(2001) in Fig. 5 and of Shemmer et al. (2000) and assuming 
an intrinsic circular symmetry of the ring, we estimate the ring 
axis inclination to the line of sight as i ~ 30°. From the appar¬ 


ent radial velocity gradient of » 20 km/s, we roughly estimate 
the deprojected physical expansion velocity as ~ 40 km/s. The 
physical dimension of the nebula is measured from the major 
axis along the N-S direction as a radius of « 30 pc (9" ). From 
these numbers we find a dynamical lifetime of the nebula of 
G ~ 8 • 10^ years. This estimate indicates a main-sequence or 
a pre-LBV status of the star. The nebula is probably a resid¬ 
ual MS bubble. Using the average density of the nebula found 
above (n^, ~ 10 cm“^), we estimate that a total mass of the neb¬ 
ula is about a few thousand solar masses. Such a big nebula 
may be composed mainly of swept up interstellar gas. 

The photometry of this star has been studied well 
(Shemmer et al. 2000), but it does not show the variability 
that is typical of LBVs. Using the Sternberg Institute plate col¬ 
lection, Zharova (2004) has found that for the last 50 years a 
variability of B416 was notably less than 1 magnitude. There 
are instead minor long-term (1987 - 1997) variations of about 
0”15 (Shemmer et al. 2000). 





Fabrika et al.: Crowded Field 3D Spectroscopy of LBV candidates 


11 


Spectrally, B416 has to be classified as an ‘iron-star” 
(Walborn and Fitzpatrick 2000). It is the brightest supergiant 
in M33 with UV-excess (Massey et al. 1996). Based on spec¬ 
tral criteria, luminosity, and the absence of a large-amplitude 
variability, B416 may be classified as B[e]-supergiant star 
(Zickgraf et al. 1986; Lamers et al., 1998b). The spectrum 
is typical of [Be]-supergiants (or ‘sg[Be]”stars, Lamers et 
al. 1998b) — very bright hydrogen lines, a forest of narrower 
Fe II lines. He I lines, forbidden lines. 

The presence of anisotropical mass loss is a common 
and distinguishing property of B[e]-supergiants (Lamers et 
al., 1998b). UV resonance line profiles in the edge-on inclina¬ 
tion stars show that the stellar winds are very slow in massive 
B[e]-supergiants (Zickgraf 2003). The expansion velocities are 
~ 70 - 100 km/s, a factor of 10 less than usually observed in 
stars of a similar spectral type. Such velocities are also typical 
of disk-like outflows in close binaries, where a donor star over¬ 
fills its critical Roche lobe. Binaries among B[e]-supergiants 
were reviewed by Zickgraf (2003), while Sholukhova et al. 
(2004) have studied radial velocities in B416 and found the 
periodical variability with a period ~ 16 days. We postpone a 
discussion of B 416 as a close binary for a future paper. 

The second star of our study, v532 shows a typical LBV- 
like variability (Romano 1978; Kurtev et al. 2001; Sholukhova 
et al. 2002). Over the past 50 years it has changed its brightness 
over a range of B = 16.3 18.4. There are large amplitude vari¬ 
ations with a time-scale of years, onto which a sporadic 0™3 - 
0™5 amplitude variability is superimposed with a time-scale 
of several months. The last irregular minimum occurred from 
1974 to 1981, and the last (also quite irregular) maximum was 
observed on 1990 - 1995. At the time of our spectroscopic ob¬ 
servations in September 1998 the star had intermediate, gradu¬ 
ally decreasing brightness. We estimate the brightness of v 532 
during this epoch as B = 17™5 +0™15 (Sholukhova et al. 2002). 

During an intermediate brightness stage in 1998 - 2001, 
the star went into a higher excitation Ofpe/WN9-like spec¬ 
trum (Fig. 10). The spectrum is closest to those of WN 10- 
11 stars, which were introduced by Smith et al. (1994) and 
studied extensively in other papers (Crowther et al. 1995abe; 
Crowther and Smith 1997; Bohannan and Crowther 1999). 
The most similar that we could identify in the literature are 
the well-known LBV AG Car during minimum (Hutsemehers 
and Kohoutek 1988, Stahl 1986, Smith et al. 1994), He 3-519 
(Smith et al. 1994), and S 142 (Crowther and Smith 1997). 
The major criteria for a very late WN 10-11 classification 
(Crowther and Smith 1997) are confirmed in the spectrum of 
v532: i) a rich low excitation nitrogen spectrum, where Nil 
lines are brighter than those of N III; ii) narrow He I lines, faint 
Hell 44686; iii) forbidden [Nil] 45755 and [Felll] lines; iv) 
the expansion velocities derived from He I P Cyg-like profiles 
are low, Voo = 100 300 km/s (it is 220 km/s in v532); v) the 

presence of a nebula with a low expansion velocity on the order 
of a few tens of km/s (marginally detected in this study). 

The spectrum of v532 is similar to the spectra of 
Ofpe/WN9 stars. These stars are believed to be directly re¬ 
lated to LBVs as in fact, “dormant” LBVs. Transitions of 
LBV <-> Ofpe/WN9 have been observed and are well-known, 
for example in AG Car, or R 127 (Stahl et al. 1983; Stahl 1986). 


We saw the same transition in v532, when during the peak of 
the maximum phase (fall 1992) the star presented a low excita¬ 
tion LBV-like spectrum (Szeifert 1996; Fabrika 2000). 

The star v 532 thus has to be classified thus as an LBV ob¬ 
ject based on its spectrum, spectral, and photometrical variabil¬ 
ity. In its current state the star may join the group of WN 11 
stars. The nitrogen spectrum of v532 indicates an evolved sta¬ 
tus, where nitrogen produced in the core appears at the surface. 

The B416 spectrum and its nebula age (~ 8-10^ years) 
indicate an MS or pre-LBV status (hydrogen-burning) of the 
star (Schaerer et al. 1993). The disk-like geometry of the neb¬ 
ula around B416 does not contradict the main sequence status 
of B 416. An asymmetry in the mass loss might be due to fast 
rotation or close binarity. 
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